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Abstract 


A search for a charged Higgs boson is performed with a data sample correspond¬ 
ing to an integrated luminosity of 19.7 ± 0.5 fb 1 collected with the CMS detector in 
proton-proton collisions at i/s = 8 TeV. The charged Higgs boson is searched for in 
top quark decays for jh h ± < m t — n%, and in the direct production pp —> t(b)H ± 
for m H ± > m t — m^. The H ± —> Uv-t and H ± —> tb decay modes in the final 
states Th+jets, }iT\ lr f+jets, and It’ (/:' =e, //) are considered in the search. No sig¬ 
nal is observed and 95% confidence level upper limits are set on the charged Higgs 
boson production. A model-independent upper limit on the product branching 
fraction B {t —> H ± b)t5(H ± —> r ± i/ T ) = 1.2-0.15% is obtained in the mass range 
7 /z h ± = 80-160 GeV, while the upper limit on the cross section times branching 
fraction u(pp —* t(b)H ± ) B{ H ± —> r ± v T ) = 0.38-0.025pb is set in the mass range 
7 jz h+ = 180-600 GeV. Here, cr(pp > Eb)H ) stands for the cross section sum 
cr(pp —> t(b)H + ) + cr(pp —> t(b)H^). Assuming B{Yi ± —> tb) = 1, an upper limit 
on cr(pp —* t(b)H ± ) of 2.0-0.13 pb is set for i h h ± = 180-600 GeV. The combination 
of all considered decay modes and final states is used to set exclusion limits in the 
77Z H ±-tan /3 parameter space in different MSSM benchmark scenarios. 
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1 Introduction 

In 2012, a neutral boson with a mass of approximately 125 GeV was discovered by the CMS and 
ATLAS experiments [THU at the CERN LHC. The properties of the new boson are consistent 
with those predicted for the standard model (SM) Higgs boson [SHU. Models with an extended 
Higgs sector are constrained by the measured mass, CP quantum numbers, and production 
rates of the new boson. The discovery of another scalar boson, neutral or charged, would 
represent unambiguous evidence for the presence of physics beyond the SM. 

Charged Higgs bosons are predicted in models including at least two Higgs doublets. The 
simplest of such models are the two-Higgs-doublet models (2HDM) IfTOl . Two Higgs doublets 
result in five physical Higgs bosons: light and heavy CP-even Higgs bosons h and H, a CP-odd 
Higgs boson A, plus two charged Higgs bosons H ± . Throughout this paper, charge conjugate 
states are implied, the cross section n(pp —> t(b)H + ) denotes the sum <r(pp —> t(b)H + ) + 
<j(pp —>• t(b)H ), and the branching fractions B( H 1 —* X) stand for B( H ± —> X). The minimal 
supersymmetric SM (MSSM) flTTHT8l used as a benchmark in this paper is a special case of a 
Type-II 2HDM scenario. In such a scenario, the couplings of the charged Higgs boson to up- 
type quarks is proportional to cot while the charged Higgs boson couplings to the down-type 
quarks and charged leptons are proportional to tan where tan [i is defined as the ratio of the 
vacuum expectation values of the two Higgs boson doublet fields. 

If the mass of the charged Higgs boson is smaller than the mass difference between the top and 
the bottom quarks, m H + < m t — n%, the top quark can decay via t —> H + b. In this case, the 
charged Higgs boson is produced most frequently via tt production. In the MSSM scenarios 
considered, it preferentially decays to a r lepton and the corresponding neutrino, H + —> r + v T , 
for tan yS > 5 Ifl9| . A representative diagram for the production and decay mode for a low-mass 
charged Higgs boson is shown in Fig. |l] (left). Compared to the SM prediction, the presence of 
the H + —> t + v t decay modes would alter the r yield in the decays of tt pairs. 

The Large Electron-Positron (LEP) collider experiments determined a model-independent lower 
limit of 78.6 GeV on the H + mass H20ti23l at a 95% confidence level (CL). The most sensitive 95% 
CL upper limits on6(t —> H + b) have been determined by the ATLAS and CMS experiments 
and are described in the following. For the H + —> r 1 v T decay mode with the hadronic decay 
of the r lepton (th) and hadronic W boson decays (ly-t-jets) final state, 95% CL upper limits 
of 1.3-0.2% have been set on £>(t —* H^b)£>(H + —> r + v T ) for m H + = 80-160 GeV by the 
ATLAS experiment using data at y/s = 8TeV |24| . For the (£=e, }i) final states 95% CL 
upper limits of 3-9% have been set by the ATLAS and CMS experiments on B(t —> H + b) in 
the H + —> r + v T decay mode for m H + = 80-160GeV assuming B( H + —> r + v T ) = 1 and using 
data at i/s = 7TeV [251126]. The H + —> cs decay mode, whose branching fraction dominates 
for tan yS < 5, has been studied by the ATLAS experiment, with 95% CL upper limits of 5-1% 
set on B( t —» H + b) for m H + = 90-160 GeV, under the assumption B( H + —> cs) = 1 and using 

data at v^s = 7TeV EZl- 

If the charged Higgs boson mass exceeds the mass difference between the top and bottom 
quark, rn H > m t — m b , the charged Higgs boson is predominantly produced by the fusion of 
bottom and top quarks illustrated in Figs. [T] (middle) and (right) for the four-flavour scheme 
(4FS) and the five-flavour scheme (5FS), respectively. In the 4FS, there are no b quarks in the 
initial state, causing a different ordering of the perturbative terms at any finite order between 
the 4FS and 5FS ll28ti3TI . The predictions of the 4FS and the 5FS cross sections calculated at 
next-to-leading order (NLO) are combined using the "Santander matching scheme" 1321 . In the 
MSSM benchmark scenarios considered, the H + —> r + v T decay mode dominates for /w H + < 
220 GeV Ifl9fl , and for large values of both m H + and tan /3, the decay H + —> tb becomes dominant 
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1 Introduction 


but the H + —> t + v t decay mode still contributes. For the H + — > r + v T decay mode, considering 
the final state with hadronic r lepton and associated W boson decays, the current upper limits 
of 0.8-0.004 pb have been set on cr (pp —> t(b)H + ) B( H + —* t + v t ) by the ATLAS experiment 
for m H + = 180-1000 GeV using data at i/s = 8TeV l2il . 



Figure 1: Left: A representative diagram for the production mode of the light charged Higgs 
boson through tt production with a subsequent decay to the Th+jets final state. Middle and 
right: Representative diagrams for the direct production of the charged Higgs boson in the 
four-flavour scheme and five-flavour scheme, respectively. 


In this paper, a search for the charged Higgs boson is performed in pp collisions at ,/s = 8 TeV. 
The data were recorded by the CMS experiment at the LHC and correspond to an integrated 
luminosity of 19.7 ± 0.5 fb 1 . The charged Higgs boson decay modes and final states discussed 
in this paper are summarized in Table [lj Model-independent limits without any assumption 
on the charged Higgs boson branching fractions are calculated on Bit —? H + b) B( H + —> r + v T ) 
and c(pp —> t(b)H + ) B( H + —> r + v T ) for m H + < m t — mt, and m H + > ni t — nit >, respectively, 
with the analysis on the H + —> r 1 v T decay mode in the Th+jets final state. Additionally, the 
H + —> t 1 v T and H + -+ tb decay modes are inclusively studied in the ftT h , single lepton (1+jets), 
and ££' {£' referring to the possible different flavour between the two leptons) final states for 
m H + > nit ~ nit,- Combined limits for the H + —> tb decay mode are set on 0 "(pp —» t(b)H + ) by 
assuming either B{ H + —>■ t + v t ) = 1 or B( H + —t tb) = 1. The T h +jets final state is not sensitive 
to the presence of charged Higgs boson decay modes other than H + —t t + v t , because any such 
decay mode would be estimated inclusively with the background through the measurement 
from data described in Section 5.2.1 All the decay modes and final states considered are used 
to set exclusion limits in the ni\\ -tan /3 parameter space for different MSSM benchmark scenar¬ 
ios ||29}[33l. To set these limits, the specific branching fractions predicted by those benchmark 
scenarios are applied. This paper includes the first results on the direct charged Higgs boson 
production for 7?z H + > rn t — nit, i n the H + —> tb decay mode. 


Table 1: Overview of the charged Higgs boson production processes, decay modes, final states, 
and mass regions analysed in this paper (£ = e, ft). All final states contain additional jets from 
the hadronization of b quarks and missing transverse energy from undetected neutrinos. The 
index after each signature denotes the section where it is discussed. 


Decay mode 

Signatures for m H + < m t — nit, 

Signatures for m H + > m t — nit, 

H+ -+ t + v t 
H+ -+ tb 

pp -+ tt -+ bH+bH'/bH+bW- 
Th+jets© 

T h+: 

firl 

pp -+ f(b)H+ 
ets©, /<Tj©, ££'0 

3, ^'0, f+jets0 


The CMS detector is briefly described in Section [2] followed by details of the event reconstruc¬ 
tion and simulation in Sections [3] and [4j respectively. The event selection together with the 
background estimation is described in Sections [5j[6j[7j and[8]for the Th+jets, }i Th, £$!, and £+jets 
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final states, respectively. The treatment of statistical and systematic uncertainties is described 
in Section [9] The results are presented in Section [To] and summarized in Section |TTj 

2 The CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diam¬ 
eter, providing a magnetic field of 3.8 T. Within the superconducting solenoid volume are a 
silicon pixel and strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and 
a brass and scintillator hadron calorimeter (HCAL), each composed of a barrel and two end- 
cap sections. Muons are measured in gas-ionization detectors embedded in the steel flux-return 
yoke outside the solenoid. Forward calorimeters extend the pseudorapidity coverage provided 
by the barrel and endcap detectors up to |//| < 5. The first level (LI) of the CMS trigger sys¬ 
tem, composed of custom hardware processors, uses information from the calorimeters and 
muon detectors to select the most interesting events in a fixed time interval of less than 4 jis. 
The high-level trigger processor farm further decreases the event rate from around 100 kHz to 
around 1 kHz, before data storage. A more detailed description of the CMS detector, together 
with a definition of the coordinate system used and the relevant kinematic variables, can be 
found in Ref. It34l . 


3 Event reconstruction 

In the data collected during 2012, an average of 21 proton-proton interactions occurred per 
LHC bunch crossing. To select the primary interaction vertex, the squared sum of the trans¬ 
verse momenta of the charged-particle tracks, £ p\, associated with each interaction vertex is 
calculated. The interaction vertex with the largest £ pj value is taken as the primary interaction 
vertex in the event l35l . The other pp collisions are referred to as pileup. 

Events are reconstructed with the particle-flow (PF) algorithm | [36l [37lj , which combines infor¬ 
mation from all sub-detectors to identify and reconstruct individual electrons, muons, photons, 
and charged and neutral hadrons. Electrons are reconstructed from clusters of ECAL energy 
deposits matched to hits in the silicon tracker It38l . Muons are reconstructed by performing a 
simultaneous global track fit to hits in the silicon tracker and the muon system lf39l . The energy 
of photons is directly obtained from the ECAL measurement, corrected for zero-suppression 
effects. The energy of charged hadrons is determined from a combination of their momentum 
measured in the tracker and the matching ECAL and HCAL energy deposits, corrected for 
zero-suppression effects and for the response function of the calorimeters to hadronic showers. 
Finally, the energy of neutral hadrons is obtained from the corresponding corrected ECAL and 
HCAL energy. The composite physics objects, such as jets, hadronic tau lepton decays, and 
missing transverse energy are reconstructed from these PF particles. 

Jets are reconstructed from the PF particles clustered by the anti-/c L algorithm [40, 41] with a 
distance parameter of 0.5. The jet momentum is determined as the vectorial sum of all particle 
momenta in the jet, and is found in the simulation to be within 5-10% of the true momentum 
over the whole pj spectrum and detector acceptance. An offset correction is applied to take 
into account the extra energy clustered in jets arising from pileup. Jet energy corrections are 
derived from simulation, and are confirmed by in situ measurements of the energy balance 
in dijet and photon+jet events H42l . Additional selection criteria are applied to each event to 
remove spurious jet-like features originating from isolated noise patterns in certain HCAL re¬ 
gions. Jets originating from pileup interactions are removed by a multivariate jet identification 
algorithm l43l . 
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3 Event reconstruction 


Jets from the hadronization of b quarks are identified (b tagged) with the "combined secondary 
vertex" algorithm [44, 45 j. The algorithm consists of evaluating a likelihood-based discrim¬ 
inator which uses information from reconstructed decay vertices of short-lived mesons and 
transverse impact parameter measurements of charged particles. In the r h +jets final state, the 
algorithm is used to identify b-tagged jets with a mistagging probability, i.e. the probability 
that a jet from the fragmentation of light quarks (u, d, s, c) or gluons is misidentified as a b 
jet, of approximately 0.1% ("tight" working point). In the analyses of the prh and f+jets final 
states, the b tagging algorithm used has a mistagging probability of 1% ("medium" working 
point), since the multijet background is smaller than in the Th+jets final state. In the analy¬ 
sis with the final state, the b tagging working point is adjusted to allow a 10% mistagging 
probability to enhance signal acceptance since the multijet background in this analysis is even 
smaller. The corresponding probability to identify a b jet is about 50, 70, and 85%, respectively. 
The difference in b tagging efficiency between data and simulation is corrected by applying 
data-to-simulation scale factors dependent on the jet px and the jet pseudorapidity (;/). 

The missing transverse momentum vector p™' ss is defined as the projection of the negative 
vector sum of the momenta of all reconstructed PF particles in an event onto the plane perpen¬ 
dicular to the beams. Its magnitude is referred to as E™ 1SS . The Ey ISS reconstruction is improved 
by propagating the jet energy corrections to it. Further filter algorithms are used to reject events 
with anomalously large E™ ISS resulting from instrumental effects Il46l . 

The "hadron-plus-strips" algorithm H47j is used to reconstruct hadronically decaying r leptons. 
The algorithm uses the constituents of the reconstructed jets to identify individual r decay 
modes with one charged and up to two neutral pions, or three charged pions. The neutral 
pions are reconstructed by clustering the reconstructed photons in narrow strips along the 
azimuthal angle direction taking into account possible broadening of calorimeter depositions 
from photon conversions. The ly candidates compatible with electrons or muons are rejected. 
Jets originating from the hadronization of quarks and gluons are suppressed by requiring that 
the r h candidate is isolated as described below. The r h identification efficiency depends on py h 
and p Th , and is on average 50% for Pj > 20GeV with a probability of approximately 1% for 
hadronic jets to be misidentified as a r h . 

Electrons, muons, and hadronically decaying r leptons are required to be isolated from other 
particles by considering transverse momenta of neutral and charged particles in a cone A R = 
(At/5) 2 + (Ap) 2 , where (p is the azimuthal angle, around the charged lepton candidate mo¬ 
mentum direction. The isolation variable for electrons, muons, and ly is defined as: 


I* = Pi + maX ( °' I] FT + I] FT - Pneutral ^eff. , 

charged \ neut. hadr. 7 / 

1,1 = I] Ft + max I 0, Pt + E] Pt - 0.5 JE p T 

charged \ neut. hadr. 7 charged,pileup / 

I Th = E FT + max ( 0/ 5E Pt “ °- 46 I] PT ) / 

charged \ 7 charged.pileup / 


( 1 ) 

( 2 ) 

(3) 


where Charged Pt is the scalar sum of the transverse momenta of charged hadrons, electrons, 
and muons originating from the primary interaction vertex, and £ neut hadr. Pt and px are 
the scalar sums over neutral hadron and photon transverse momenta, respectively, in the cone 
A R around the charged lepton candidate momentum direction. The presence of particles from 
pileup events is taken into account depending on the charged-lepton type. For electron can¬ 
didates, the scalar sum of the px of photons and neutral hadrons from pileup events in the 
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isolation cone is estimated as the product of the neutral-particle transverse momentum density 
and the effective cone area, p n eutrai /left.- The p n eutrai component is evaluated from all photons 
and neutral hadrons in the event, and A e ff accounts for the presence of pileup events. For 
muons and hadronically decaying t leptons, the scalar sum of the pj of photons and neutral 
hadrons from pileup events is estimated from the scalar sum of the transverse momenta of 
charged hadrons from pileup events in the isolation cone, ^charged,pileup Ptv by multiplying it 
by the average ratio of neutral- to charged-hadron production in inelastic pp collisions. Since 
the contribution from neutral hadrons is ignored when computing the Ty isolation variable, 
the pileup correction factor is slightly smaller than that used for correcting the muon isolation 
variable. 

For electrons, an isolation cone size of A R = 0.3 or 0.4 is used, depending on the final state. 
For muons and hadronically decaying r leptons, isolation cone sizes of A R = 0.4 and 0.5 are 
used, respectively. Electrons and muons are considered isolated if the relative isolation variable 
7^1 = 7*7 p^, where i = e, p, is lower than 10-20%, depending on the final state. Hadronically 
decaying r leptons are considered isolated if I Th < 1 GeV. 


4 Simulation 

The signal processes are generated with pythia 6.426 [48| . The tt, W+jets, and Z+jets back¬ 
grounds are generated using the MadGraph 5.1.3.30 l49l event generator with matrix ele¬ 
ments (ME) providing up to four additional partons, including b quarks. The event generator 
is interfaced with PYTHIA to provide the parton showering and to perform the matching of the 
soft radiation with the contributions from the ME. The single top quark production is generated 
with POWHEG 1.0 Il50l - l54i and the quantum chromodynamics (QCD) multijet and diboson pro¬ 
duction processes WW, WZ, and ZZ are generated using pythia. Both the MadGraph and 
POWHEG generators are interfaced with PYTHIA for parton shower and hadronization. The 
TAUOLA 27.121.5 l55l package is used to generate r decays for the simulated signal, as well as 
background samples. 

The events are passed through full CMS detector simulation based on Geant4 [56] 157], fol¬ 
lowed by a detailed trigger simulation and event reconstruction. Simulated minimum bias 
events are superimposed upon the hard interactions to match the pileup distribution observed 
in data. The PYTHIA parameters for the underlying event are set according to the Z2* tune, 
which is derived from the Z1 tune I f58l , which uses the CTEQ5L parton distribution set, whereas 
Z2* adopts CTEQ 6 L E3. 

The number of tt events produced is normalized to the predicted tt production cross section of 
246.7lg4 ± 11.4 pb as calculated with the TOP++ v2.0 program to next-to-next-to-leading order 
(NNLO) in perturbative QCD, including soft-gluon resummation to next-to-next-to-leading- 
logarithmic (NNLL) order [60[ , and assuming m t = 173.34 GeV [MT| . The first uncertainty orig¬ 
inates from the independent variation of the factorization and renormalization scales, \i and 
Pr, while the second is associated with variations in the parton density functions (PDFs) and 
strong coupling constant as, following the PDF4LHC prescription with the MSTW2008 68 % 
CL NNLO, CT10 NNLO and NNPDF2.3 5-flavour fixed-flavour number (FFN) PDF sets [62] - 
[65| . The predicted cross section is in good agreement with the measurements by ATLAS and 
CMS l66l[67l . The top quark pj spectrum in data is found to be softer than that predicted using 
the MadGraph MC generator l 68 l . To correct for this effect, the tt events are reweighted to 
make the top quark pj spectrum in simulation match that observed in data [691 . 

The NNLO SM prediction is calculated with FEWZ v3.1 for the W+jets and Z/ 7 * backgrounds [70l 
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5 The T h +jets final state for H+ -» r+v T 


IZD- The cross section for the t-channel single top quark sample is calculated at next-to-leading 
order (NLO) in QCD with HATHOR v2.1 Il72l l73il with PDF and ttg uncertainties calculated using 
the PDF4LHC prescription (62] [741 . For the single top quark s-channel and tW-channel cross 
section, the SM prediction at NNLL in QCD is taken from Refs. 17511761. 


5 The Th+jets final state for H+ — ► r + v T 

In this analysis, a charged Higgs boson is assumed to be produced through the tt —> bH + bH“, 
tt —> bH + bW% and pp —> t(b)H + processes and searched for in the H + — > t + v t decay mode 
with a hadronic decay of the r and a hadronic decay of the W boson that originates from the 
associated t —? bW decay In these events, the missing transverse momentum is expected 
to originate from the neutrinos in the decay of the charged Higgs boson, which allows the 
reconstruction of the transverse mass, nij, of the charged Higgs boson: 


nij = iss (l - cos A (p(pj h , p™ iss )), (4) 

where Py 1 denotes the transverse momentum vector of the hadronically decaying r lepton and 
Pj' its magnitude, and Atp is the angle between the ip direction and the p| TI,ss in the transverse 
plane. The presence of the two neutrinos from the charged Higgs boson decay smears the 
expected Jacobian peak somewhat, but leaves the kinematic edge at the charged Higgs boson 
mass intact. The search is performed as a shape analysis, using the transverse mass to infer the 
presence of a signal. The dominant background processes are the SM tt and single top quark 
production, and the electroweak (EW) processes: W+jets, Z+jets, and dibosons (WW, WZ, ZZ). 
The multijet background constitutes a subleading background. 

5.1 Event selection 

Events are selected with a trigger that requires the presence of a q, and large Ey ,ss . First the 
events are required to have calorimetric Ey ISS > 40GeV at the first level of the CMS trigger 
system. The calorimetric E™ 1SS is defined as the Ey ISS calculated from the ECAL and HCAL 
energy deposits instead of the PF particles. At the high-level trigger, the events are required to 
have calorimetric E™ 1SS > 70GeV, and a r h of py > 35GeV and ;/ Th | < 2.5. The q, is required 
to be loosely isolated, to contain at least one track of pj > 20 GeV, and to have at most two 
tracks in total, targeting the r lepton decays into a single charged pion and up to two neutral 
pions. The probability for a signal event to be accepted by the trigger amounts to 8-14% in the 
m H + range of 80-160 GeV, and 19-44% in the »z H range of 180-600 GeV with all tau decays 
considered. 

The efficiency of the r part of the trigger is evaluated as a function of py 1 ' using a "tag-and- 
probe" technique It47l from Z/ 7 * -+ q,T h events, where q, refers to a muonic r lepton decay. 
The efficiency of the Ey ISS part of the trigger is evaluated from events with a tt-like final state of 
Th+jets selected with a single-r trigger. The trigger efficiencies in simulated events are corrected 
with data-to-simulation scale factors applied as function of Py for the r h part of the trigger and 
as function of Ey ISS for the Ey ISS part of the trigger. The scale factors range between 0.95-1.06 
and 0.97-1.02 for the Th and Ey ISS parts of the trigger, respectively. 

Selected events are required to have at least one q, with py h > 41 GeV within |?/| < 2.1 and 
to be matched to a trigger-level q, object. These thresholds are chosen to be compatible with 
the single-muon trigger used for estimate of backgrounds with hadronic r decays from control 
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samples in data, as described in Section 5.2.1| Only one charged hadron is allowed to be associ¬ 
ated with the Th and its pj is required to fulfil pj > 20 GeV. Background events with W —> rv T 
decays are suppressed by requiring R T = ^charged hadrony pT h > g j ypg observable is sen¬ 
sitive to different polarizations of r leptons originating from decays of W bosons (spin 1) and 
from decays of H + (spin 0) If77| . 


A tt-like event topology is selected by requiring at least three jets of pj > 30 GeV and ;/| < 
2.4 in addition to the r h and by requiring at least one of the selected jets to be identified as 
originating from the hadronization of a b quark. To select a fully hadronic final state, events 
containing identified and isolated electrons (muons) with pj > 15 (10) GeV are rejected. The 
electron (muon) candidates are considered to be isolated if the relative isolation E® el (/ ”, ( ), as 
described in Section[3} is smaller than 15% (20%). 

To suppress the multijet background, Ey ISS > 60 GeV is required. The lower Ey ISS threshold on 
the PF Ej' ss compared to the calorimetric E“ 1SS requirement applied at the high-level trigger 
improves the signal acceptance for m H + < m t — n%. This approach can be used because the PF 
E™ ISS has better resolution than the calorimetric Ey ISS [|46| . 

In the multijet events selected with the r+Ey ISS trigger a hadronic jet is misidentified as the 
Th in the event. In addition, the Ty typically has a recoiling jet in the opposite direction. The 
Ey ISS in these events arises from the mismeasurement of the momenta of these jets with the 
^miss di rec tion aligned with py. The best performance for multijet background suppression 
and signal acceptance is obtained with two-dimensional circular selections instead of simple 
selections based on azimuthal angle differences. The variables used for the azimuthal angle 
selections are defined as 


R rou = min | yj ( A( H T h/ PT iss )Y + {n~ A<^(jet n ,p^ iss )) 2 |, 
R£T = min | sj{n - A^(r h , p™*)) 2 + (A^(jet„, pf ss )) 2 } , 


(5) 


where the index n refers to any of the three highest pj jets in the event and A cp denotes the 
azimuthal angle between the reconstructed p™' ss and the ip or one of the three highest-px jets. 
The labels "coll" and "bb" denote the collinear and back-to-back systems of the Ty and the Ey ISS , 
respectively The selected events are required to satisfy > 0.70 and Kj^ 1 ' 1 > 0.70. 

The same event selection is used for all the values considered. 


5.2 Background measurements 

The background contributions arise from three sources: 


1. Irreducible background from EW processes — W+jets, Z+jets, and dibosons — as well as 
SM tt and single top quark production, where the selected t h originates from a hadronic 
decay of a r lepton ("EW+tt with Th"). 

2. Reducible background from multijet events with large mismeasured Ey ISS and jets that 
mimic hadronic r decays. 

3. Reducible background from EW+tt events, where an electron, muon, or a jet is misiden¬ 
tified as the Th ("EW+tt no Th"). 
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5 The T h +jets final state for H+ -» r+v T 


The two largest backgrounds, "EW+tt with Th" and multijets, are measured from control sam¬ 
ples in data, as explained in Sections 5.2.1 and 5.2.2 The contribution from "EW+tt no Th" is 
estimated from simulation and is described in Section l5.2.3l 


5.2.1 Measurement of the EW+tt with hadronically decaying r leptons background 

The nij distribution for the "EW+tt with Th" background is modelled via an embedding tech¬ 
nique. It uses a control data sample of p+jets events selected with a single-// trigger. The same 
jet selection as in the Th+jets sample is used, and events with electrons or additional muons 
are rejected. Then, the selected // is replaced by a simulated r lepton decay. The simulated t 
lepton momentum is the same as that of the selected //, and the reconstructed t decay products 
are merged with the original //+jets event, from which the reconstructed muon is removed. In 
these hybrid events, the jets are reclustered and the Ey ISS is recalculated and then the events 
are subjected to the same event selection as the Th+jets sample, i.e. Th identification, b tagging, 
Ey ISS requirement, and the azimuthal angle selections are applied. 

To obtain the nij distribution for the "EW+tt with Th" background, the effect of the muon trigger 
and the muon offline reconstruction need to be unfolded, and the efficiency of the T+Ey ISS 
trigger must be taken into account. First, the weight of each hybrid event is increased by the 
inverse of the muon trigger and identification efficiencies. Then, the efficiency of the T+Ey ISS 
trigger is applied by weighting the events with the efficiencies of the t part of the trigger and 
the first trigger level part of the E™ ISS trigger. The rest of the Ey ,ss part of the trigger is taken 
into account by applying a requirement on a hybrid calorimetric Ey ISS constructed from the 
original event and the simulated t lepton decay. 

After the trigger has been taken into account, further corrections are applied. In a fraction of 
the selected //+jets events the // originates from a decay of a t lepton, leading to an overestima¬ 
tion of the EW+tt background by a few percent. This bias is corrected for by applying to the 
hybrid events py-dependent correction factors derived from simulated tt events. A residual 
difference is seen in the nij distribution between non-embedded T+jets and embedded //+jets 
events in simulated tt events. This difference is corrected by weighting the hybrid events by 
mj-dependent correction factors derived from simulated tt events. The tt events constitute 
about 85% of the "EW+tt with T h " background. 

It should be noted that the embedding technique allows the separation of signal from the H ! -+ 
t + v t decay mode from other decay modes, such as H + -+ tb, where the t lepton originates 
from a W boson decay. Namely, in the other charged Higgs boson decays, x leptons and muons 
are produced at equal rates causing the embedding technique to include the H + —> tb signal 
from data (and other such signals) as part of the "EW+tt with Th" background. 


5.2.2 Measurement of the multijet background 


The multijet background is measured with a "Th misidentification rate" technique. An estimate 
of the multijet background is obtained by measuring the probability of the Th candidate to pass 
the nominal and inverted Th isolation criterion. The misidentification rate is measured in bins of 
Th transverse momentum, in an event sample that is obtained pri or to applying the b tagging, 
Ey ISS , and parts of the event selection described in Section 5.1 The event sample that 


passes the nominal Th isolation selection contains a nonnegligible contamination from EW+tt 
backgrounds with genuine and misidentified t leptons. Therefore, the number of multijet and 
EW+tt events is determined by a maximum likelihood fit of the Ey ISS distribution. A fit is 
performed for each py h bin. For multijet events, the Ey ISS templates are obtained from the 
data sample with inverted Th isolation by subtracting a small contribution of simulated EW+tt 
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events. The E™ ISS templates for the EW+tt events are taken from simulation in the nominal 
region. The misidentification rate probabilities Wj are defined as the ratio of the number of 
multijet events in the isolated sample and the inverted isolation sample. Their measured values 
vary between 0.050-0.061 depending on the Pj bin with a statistical uncertainty smaller than 

3%. 

The measured T b misidentification rate probabilities are then applied as weights to multijet 
events passing all nominal event selection criteria, except that the isolation criterion applied 
on the Ty is inverted. The number of multijet events is obtained by subtracting the number of 
simulated EW+tt events from data. The estimate for the number of multijet events in a given 
bin i of the vij distribution (N mLlltl|lt ) is obtained by summing these weighted events over the 
Pj' bins according to 


,rmultijet \ 1 / a rdata. inverted a r EW+tt, inverted ir\ 

N i =h( N i,j ~ N i,j > w i’ ( 6 ) 

i 

where N is the number of events and i and j denote nij and Pj' bins, respectively 

5.2.3 The EW+tt with misidentified r leptons background 

The "EW+tt no Th" background originates almost solely from jets that are misidentified as the 7}, 
with a small contribution from electrons and muons misidentified as the Ty. About 85% of the 
"EW+tt no r h " background events come from tt and the rest from single top quark production 
in the tW- and t-channels. The number of selected simulated events in the single top quark 
samples is small and therefore the nij distribution for them is estimated with a procedure where 
the probability of each event to pass the b tagging is applied as a per-event weight instead of 
applying the b tagging selection. This probability is evaluated for simulated events with the 
tt-like final state as function of jet pj and flavour. 

5.3 Event yields 

Figure [2] shows the event yields after each selection step starting from the requirement that a 
Th, no isolated electrons or muons, and at least three jets are present in the event. The multijet 
background and the "EW+tt with Th" background are shown as measured from data, while the 
"EW+tt no Th" background is shown as estimated from the simulation. The data agree with the 
sum of expected backgrounds within the total uncertainties. 

The observed numbers of events after the full event selection are listed in Table |2j along with 
those expected for the backgrounds and for the charged Higgs boson production. The sys¬ 
tematic uncertainties listed in Tableware discussed in Section [9] The nij distributions with all 
event selection criteria applied are shown in Fig. [3] for m H + < m t — and m H + > ;» t — v%. 
In the 77Z h + > rn L — m b region, the limited number of background events in the high-mj tail 
is modelled by fitting an exponential function of the form poe~ p “ mT_c \ where po and p\ are 
positive free parameters and where c = 180 GeV is the starting point of the fit. In the region of 
mj > 160 GeV the event yields for the backgrounds are replaced by those obtained from this 
fit. The slight excess of observed events in the nij spectrum for m H ■ > 7/z t — and limits on 
the production of the charged Higgs boson extracted from these distributions are discussed in 
Section [Tol 
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Figure 2: The event yield in the Ty+jets final state after each selection step. For illustrative 
purposes, the expected signal yields are shown for m H + = 120 GeV normalized to 13 (t —> 
H + b) H + —> r + i/ T ) = 0.01 and for m H + = 300 GeV normalized to u(pp —» t(b)H + ) B{ H + —> 
t 1 v T ) = 1 pb, which are typical values for the sensitivity of this analysis. The bottom panel 
shows the ratio of data over the sum of expected backgrounds and its uncertainties. The cross- 
hatched (light grey) area in the upper (lower) part of the figure represents the statistical un¬ 
certainty, while the collinear-hatched (dark grey) area gives the total uncertainty in the back¬ 
ground expectation. 


6 The fir h final state for H+ -* r + t/ T and H+ —► tb 

In this analysis, a charged Higgs boson with m H + > »z t — is assumed to be produced 
through pp —> t(b)H + : this can result in a final state characterized by the presence of two 
leptons. Here we describe the //Th choice, whereas the it' (t = e, }i) final state is discussed in 
Section[7] The }lt h final state is sensitive to the charged Higgs boson decay modes H + —* t + v t 
and H + —> tb. 

In the first case, the r decays hadronically and the final state is characterized by the leptonic 
decay of the W boson from the t —> bW decay which results in a muon in the final state. In 
the second, at least one of the W bosons from the top quarks decays to a r lepton which in turn 
decays to hadrons, whereas the other decays into a muon. Selecting the tau decay for one of 
the W bosons enhances the sensitivity to the H + —> r + v T decay mode of the charged Higgs 
boson. In this final state, the charged Higgs boson production is characterized by a number of 
b-tagged jets larger than in the SM backgrounds, and consequently the shape of the b-tagged jet 
multiplicity distribution is used to infer the presence of a signal. The dominant SM background 
processes are from tt —> /zip + X, and other backgrounds where a jet is misidentified as a Th 
(mainly lepton+jet tt events and W+jet production). 

6.1 Event selection 

The event selection is similar to that used in the measurement of the top quark pair production 
cross section in dilepton final states containing a Th [78, Z?j. A single-muon trigger with a 
threshold of pj > 24GeV and \rj\ < 2.1 is used to select the events. 
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Table 2: Numbers of expected signal and background events with their statistical and system¬ 
atic uncertainties listed together with the number of observed events after the full event selec¬ 
tion is applied in the r h +jets final state. For illustrative purposes, the expected signal yields 
are shown for m H + = 120GeV normalized to B {t —> H + b)S(H + —* r + i/ T ) = 0.01 and for 
m H + = 300GeV normalized to cr (pp —* t(b)H + )S(H + —> t + v t ) = lpb, which are typical 
values for the sensitivity of this analysis. 


Source 

Nevents( ± Stat ± syst) 

Signal, 7/z h + = 120 GeV 

151 ±4+^ 

Signal, m H + = 300 GeV 

168 ±2 ± 16 

EW+tt with Th (data) 

283 ± 12 

Multijet background (data) 

80 ± 3 

EW+tt no Th (sim.) 

47 ± 2 

Total expected 

410 ± 12 

Data 

392 
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Figure 3: The transverse mass (nij) distributions in the r^+jets final state for the H + mass 
hypotheses of 80-160 GeV (left) and 180-600 GeV (right). The event selection is the same in 
both left and right plots, but in the right plot the background expectation is replaced for nij > 
160 GeV by a fit of the falling part of the nij distribution. Since a variable bin width is used 
in the left plot the event yield in each bin has been divided by the bin width. For illustrative 
purposes, the expected signal yields are shown in the left plot for 7/z H + = 120 GeV normalized 
to B(t —> H + b) B( H + —* r + v T ) = 0.01 and in the right plot for m H + = 300GeV normalized 
to f(pp —> t(b)H + ) B( H + —> r + v T ) = 1 pb, which are typical values for the sensitivity of this 
analysis. The bottom panel shows the ratio of data over the sum of expected backgrounds 
along with the uncertainties. The cross-hatched (light grey) area in the upper (lower) part of 
the figure represents the statistical uncertainty, while the collinear-hatched (dark grey) area 
gives the total uncertainty in the background expectation. 


Events are selected by requiring one isolated muon with pj > 30GeV and \rj\ < 2.1, one 
hadronically decaying r with pj > 20GeV and \rj\ < 2.4, at least two jets with pj > 30GeV 
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6 The /<r h final state for H+ -» r + v T and H+ -> tb 


and | rj \ < 2.4, with at least one jet identified as originating from the hadronization of a b quark, 
and Ej' ss > 40GeV. The q, and the muon are required to have opposite electric charges. The 
muon candidate is considered to be isolated if the relative isolation, as defined in Section [3} is 
/ re l < 0.12. The muon and the r are required to be separated from each other and from any 
selected jet by a distance A R > 0.4. The choice of the radius matches the lepton isolation cone. 
Events with an additional electron (muon) with I re \ < 0.2 and pj > 15(10) GeV are rejected. 


6.2 Background estimate 


There are three main background categories. The first includes backgrounds that contain a 
genuine muon and a genuine Tp, and is constituted by tt —> MT h + X production, associated 
tW —> px h + X production, Z -+ tt —* px h Drell-Yan production, and W —> //r h + X pro¬ 
cesses. The second category includes backgrounds with a genuine muon and an electron or 
muon misidentified as a Ty, namely tf —> pi + X, Z —» up, associated tW -+ pi + X pro¬ 
duction and VV —> pi + X production. The third category involves processes with a genuine 
muon and a jet misidentified as a Th, which include tt —> p+ jets, V+jets, single top quark, 
and VV —> »+jets events. Within those categories, all genuine muons come from W/Z de¬ 
cays, either direct (W —» pv, Z —> pp) or via intermediate r decays (W —> xv —t p + E™' ss ; 
Z —> tt —y pT ]h + E™ 1SS ). 


The backgrounds from the first two categories are estimated using simulation, except for the 
background due to Z/ 7 * —> xx events with one Ty and one r decaying into a muon, which is 
estimated by taking for each variable the normalization from simulation and the shape from 
Z -+ pp events in data, where each muon has been replaced with reconstructed particles from 


a simulated r lepton decay. The procedure is similar to the one described in Section 5.2.1 


The backgrounds containing a jet misidentified as a Ty come mostly from W+jets and from 
tf —> W 1 W bb —> ^vqq'bb events, and are collectively labeled "misidentified r h " in the fol¬ 
lowing tables and plots. This background is estimated by weighting each event in a p+ > 3 jets 
control sample by the probability for any jet in the event to mimic a Ty. The contribution from 
tf —y pi + X events, where one jet fakes a Th, is estimated using simulation and is subtracted 
from the data driven estimate to avoid double counting. The probability that a jet is misiden¬ 
tified as a Th is measured from data as a function of jet pj, rj, and jet radius using W+jets and 
multijet events [47j [79). Here, the estimate of the misidentified Th background is improved 
with respect to the method used in Ref. H25l by weighting according to the quark and gluon jet 
compositions (from simulation) the estimates obtained in the W+jet and multijet samples 11791 . 
This data driven estimate is different from the one described in Section 5.2.2[ where the control 
region is obtained by inverting isolation requirements on the reconstructed Th and only one 
control region is used. Here, estimating the fake rate in multijet events is not enough: the con¬ 
tamination from W+jets and tf —» W 1 W bb -+ pvqq'bb events must be taken into account as 
well. The improvement in the central value of the estimate is verified with a closure test con¬ 
sisting in applying the data driven method to simulated events: the result of the closure test is 
compatible with the yields obtained from simulation, within the uncertainties. The systematic 
uncertainty associated to the data driven method is reduced by 30% with respect to the cited 
paper. The misidentified Th background measured from data is consistent with the expectations 
from simulation. 


The fraction of events from SM tf production that is not included in the tf —> px h + X or 
misidentified Th contributions is labeled as "other tf" in the following tables and plots. The 
tf events are categorized in order to separate the contribution from each decay mode, using the 
full information on the simulated particles. 
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The single lepton trigger efficiency and the muon isolation and identification efficiencies are 
corrected by multiplicative data-to-simulation scale factors that depend on the muon pj and 
rj. Those factors are derived using a "tag-and-probe" method [80l[8ll . The trigger correction 
factors vary between 0.96 and 0.99, whereas the corrections to isolation and identification effi¬ 
ciency vary between 0.97 and 0.99. 

6.3 Event yields 

The numbers of expected events for the SM backgrounds, the expected number of signal events 
from the pp —> t(b)H + process for m H + = 250 GeV for the decay modes H + —> tb and H + —> 
T + v Zr and the number of observed events after all the selection requirements are summarized 
in Table [3] Statistical and systematic uncertainties evaluated as described in Section [9] are also 
shown. For illustrative purposes, the number of signal events is normalized, assuming a 100% 
branching fraction for each decay mode, to a cross section of 1 pb, which is typical of the cross 
section sensitivity of this analysis. 

Table 3: Numbers of expected events in the f/ry final state for the SM backgrounds and in the 
presence of a signal from H + — > tb and H + —> t + v t decays for m H + = 250 GeV are shown 
together with the number of observed events after the final event selection. For illustrative 
purposes, the number of signal events is normalized, assuming a 100% branching fraction for 
each decay mode, to a cross section of 1 pb, which is typical of the cross section sensitivity of 
this analysis. 


Source 

N ev ents ( ± Stat ± syst) 

H+ -> t + v T/ m H , = 250GeV 

176 ± 10 ± 13 

H+ tb, m w =250 GeV 

37 ± 2 ± 3 

tt — y jiT\^ + X 

2913 ± 14 ± 242 

Misidentified r h 

1544 ± 14 ± 175 

tt dilepton 

101 ± 10 ± 27 

Z/y* —> ee,pp 

12 ±3 ±4 

Z/y* —>■ TT 

162 ± 40 ± 162 

Single top quark 

150 ± 12 ± 18 

Dibosons 

20 ± 3 ± 2 

Total SM backgrounds 

4903 ± 45 ± 341 

Data 

4839 


Data and simulated event yields at various steps of the event selection are shown in Fig. [4] 
(left). Since the background estimate is derived from data only after requiring one r h , the 
backgrounds here are normalized to the SM prediction obtained from the simulation. A good 
agreement (~ 1% after the full selection) is found between data and the SM background ex¬ 
pectations. The multijet background contribution is negligible at the final selection step. The 
expected signal event yields are shown as dashed lines. 


The b-tagged jet multiplicity after the full event selection is shown in Fig . |4| ( right). Here the 
misidentified Th background is derived from data, as discussed in Section 6.2 The ratio of the 
data to the sum of the expected SM background contributions is shown in the bottom panel. 
Limits on the production of the charged Higgs boson are extracted by exploiting this distribu¬ 
tion. 










14 


6 The /<r h final state for H+ -» r + v T and H+ -> tb 
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Figure 4: Left: event yields after each selection step, where OS indicates the requirement to 
have opposite electric charges for the ry and the }i. The backgrounds are estimated from sim¬ 
ulation and normalized to the SM prediction. Right: the b-tagged jet multiplicity distribution 
after the full event selection. As opposed to the left plot, the "misidentified Th" component 
is estimated using the data-driven method and labeled misID (DD)", while the remaining 
background contributions are from simulation normalized to the SM predicted values. For both 
distributions, the expected event yield in the presence of the H + —y tb and H + —> r + v T decays 
is shown as dashed lines for «z H + = 250 GeV. For illustrative purposes, the number of signal 
events is normalized, assuming a 100% branching fraction for each decay mode, to a cross sec¬ 
tion of 1 pb, which is typical of the cross section sensitivity of this analysis. B( H + —> tb) = 1 
and B{ H + —> t + v t ) = 1, respectively. The bottom panel shows the ratio of data over the sum 
of the SM backgrounds; the shaded grey area shows the statistical and systematic uncertainties 
added in quadrature. 
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7 The dilepton (ee/e///////) final states for H + —» r + v T and H+ —► tb 

In this analysis, a charged Higgs boson with m H + > m t — m b is assumed to be produced 
through pp —> t(b)H + and is searched for in the ££! final state. Assuming that the top quark 
produced in association with the charged Higgs boson decays as t —>■ £vb, the dilepton final 
state is sensitive to charged Higgs boson decay modes H + —y tb (via leptonic decays of the top) 
or H + — y t~v t (via leptonic decays of the tau lepton). 

This leads to a final state similar to the SM tt dilepton final state, with the addition of one or 
two b jets. The shape of the b-tagged jet multiplicity distribution is used to infer the presence 
of a charged Higgs boson signal. The dominant SM backgrounds are from tt and single top 
quark production. An optimization procedure selected the b-tagged jet multiplicity variable as 
the most discriminating between the signal and the main backgrounds. 

7.1 Event selection 

The event selection is similar to that used for the measurement of the SM tt cross section and 
of the ratio B{ t — y Wb) / B(t —y Wq) in the dilepton channel |EE!§21- Data were collected 
with double-lepton triggers (ee/ pp / ep) with pj thresholds of 17GeV for the leading lepton 
and 8GeV for the other. After offline reconstruction, events are required to have two iso¬ 
lated, oppositely charged, leptons (one electron and one muon, or two electrons, or two muons) 
with pj > 20GeV and \t]\ < 2.5 (\p\ < 2.4) for electrons (muons), and at least two jets with 
p T > 30GeV and \rj\ < 2.4. The relative isolation requirement is l re i < 0.15(0.20) for electrons 
(muons). Jets are required to be separated by a distance A R = 0.4 from the isolated leptons. 
A minimum dilepton invariant mass of 12 GeV is required to reject SM background from low- 
mass resonances. For the same flavour channels (ee, pp), events with dilepton invariant mass 
within 15 GeV from the Z boson mass are vetoed. In order to account for the presence of neu¬ 
trinos, Ej ,ss > 40 GeV is required. Finally, at least two b-tagged jets are required. 

7.2 Background estimate 

The main background comes from tt events in which both W bosons decay leptonically, and 
surpasses by more than one order of magnitude the sum of the remaining backgrounds. All 
backgrounds are estimated from simulation. The dilepton trigger efficiency is corrected by 
a multiplicative data-to-simulation scale factor dependent on the final state, in order to pro¬ 
vide agreement between data and simulation; the corresponding scale factors are computed 
using the "tag-and-probe" method, and the resulting values are 0.97, 0.95, and 0.92 for the ee, 
ep, and pp final states, respectively. The data-to-simulation scale factors for the lepton iden¬ 
tification and isolation efficiencies are defined using a second "tag-and-probe" method with 
Z — y e'e /p + p~ events. For electrons (muons) with pj > 20GeV, they are found to vary 
between 0.91 (0.97) and 1.0 (0.99). 

7.3 Event yields 

The number of data events after each selection requirement are in good agreement with the SM 
background expectations, and are shown in Fig. [5] (left), for the ep final state as a representative 
example. 

The number of expected events after all selections in the ££’ final state is summarized in Ta¬ 
ble [4] for the SM background processes and for a charged Higgs boson with a mass of »z H + = 
250 GeV. The main background comes from tt production in the dilepton final state, includ¬ 
ing all three lepton flavours. Backgrounds from tt production in the final states other than 
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8 The single-lepton (e/p+jets) final states for H+ -» tb 
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Figure 5: The event yields at different selection steps (left) and the b-tagged jet multiplicity 
after the full event selection for the ep final state (right). For illustrative purposes, the number 
of signal events is normalized, assuming a 100% branching fraction for each decay mode, to a 
cross section of 1 pb, which is typical of the cross section sensitivity of this analysis. The bottom 
panel shows the ratio of data over the sum of the SM backgrounds; the shaded area shows the 
statistical and systematic uncertainties added in quadrature. 


"tt dilepton" (labelled "other tt") and other SM processes result in significantly smaller yields. 
Statistical and systematic uncertainties evaluated as described in Section[9]are also shown. The 
data agree with the sum of expected backgrounds within the total uncertainties. 

The b-tagged jet multiplicity distribution for the ep final state, shown after the full event selec¬ 
tion in Fig. [5] (right), is used to extract limits on the charged Higgs boson production. 


8 The single-lepton (e///+jets) final states for H+ —► tb 

In this analysis, a charged Higgs boson with m H + > m { — m b and produced in association 
with a top quark pp —> t(b)H + , is searched for in the decay mode H + —> tb. Of the two W 
bosons produced from the top quark decays, one decays leptonically, while the other decays 
hadronically, leading to the final state signature of one lepton, jets, and E™ ISS . These final states 
are similar to the SM tt semileptonic final states, with the addition of one or two b jets. While the 
dilepton analysis (Section [7j uses the shape of the full b tagged jet multiplicity distribution to 
check for the presence of a signal, for this analysis, an optimization procedure led to use of the 
Hj distribution, defined as the scalar sum of the px of all selected jets, subdivided by b tagged 
jet multiplicity, to infer the presence of a charged Higgs signal. Due to the jet composition 
of the signal, the Hj distribution peaks at higher energies and has a less steeply falling high 
energy tail than the major backgrounds. The dominant backgrounds are tt, W+jets, and single 
top quark production. 

8.1 Event selection 

Data were collected by the single-electron or a single-muon trigger with px thresholds of 27 and 
24 GeV, respectively. The offline event selection requires the presence of exactly one isolated 
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Table 4: Number of expected events for the SM backgrounds and for signal events with a 
charged Higgs boson mass of m H = 250 GeV in the ee, ep, and pp dilepton final states after 
the final event selection. For illustrative purposes, the number of signal events is normalized, 
assuming a 100% branching fraction for each decay mode, to a cross section of 1 pb, which 
is typical of the cross section sensitivity of this analysis. Event yields are corrected with the 
trigger and selection efficiencies. Statistical and systematic uncertainties are shown. 


Source 

ee 

ep 

m 

H+ -+ r+v T , m H + = 250 GeV 

39 ± 3 ± 3 

97 ± 4 ± 5 

40 ± 3 ± 3 

H+ -+ tb, m w = 250 GeV 

85 ± 3 ± 2 

219 ±5 ±5 

90 ± 3 ± 2 

tt dilepton 

5692 ± 17 ± 520 

15296 ± 28 ± 1364 

6332 ± 18 ± 572 

Other tt 

22 ± 4 ± 5 

40 ± 5 ± 9 

17 ± 3 ± 5 

Z/ 7 * -+ U 

96 ± 7 ± 35 

36 ± 2 ± 7 

139 ± 10 ± 42 

W+jets, multijets 

6 ± 2 ± 1 

3±1 ±1 

< 1 

Single top quark 

199 ±10 ±21 

522 ± 15 ± 54 

228 ± 10 ± 26 

Dibosons 

15 ±1 ±2 

43 ± 2 ± 6 

20 ± 1 ± 3 

Total SM backgrounds 

6032 ± 20 ± 521 

15941 ± 32 ± 1365 

6736 ± 23 ± 575 

Data 

6162 

15902 

6955 


electron (muon) with pj > 30 (27) GeV and \rj\ <2.5 (2.4). The electrons (muons) are required 
to be isolated with 1^ < 0.10 (0.20), with F cj defined in Section [ 3 ] Events with additional 
leptons are rejected. To maintain exclusivity with the other analyses included in this paper, 
events with one or more hadronic r decays with pj > 20GeV and \t] Th \ < 2.4 are rejected. In 
addition, the presence of at least two jets with pj > 30 GeV and |//| < 2.4 are required, with 
px > 50 GeV for the jet with the highest pj. At least one of the selected jets is required to be 
b-tagged. The E™ 1SS must exceed 20 GeV to mimic the presence of a neutrino in the final event 
signature. 

To account for differences in modelling of the lepton identification and trigger efficiency be¬ 
tween simulation and data, rj- and px-dependent scale factors are applied. The single-electron 
trigger correction factor is 0.973 (1.020) for \rj\ < 1.5 (1.5 < \t]\ < 2.5) and the single-muon 
trigger correction factors and corrections to identification efficiency are similar to those in Sec¬ 
tion \62\ 

Events are classified into two categories, a signal region (SR) and a control region (CR). The 
CR is defined by having low reconstructed jet multiplicity, 2 < Nj e t < 3, and is used to de¬ 
rive normalizations for dominant backgrounds from data. The SR is distinguished by its high 
jet multiplicity, and defined by the requirement Nj e t > 4. These categories are further subdi¬ 
vided according to the b-tagged jet multiplicities, N b ta g/ with the CR split into 3 subcategories 
(/Vb tag = 0, Nb tag = 1/ and Nb tag > 2) and the SR split into two (Nb tag = 1 and Nb ta g > 2). 
Distinguishing between electron and muon channels leads to a total of four SR categories and 
six CR categories. 

8.2 Background estimate 

The following background processes are considered: tt, W+jets, single top quark, Z/ 7 *+jets, 
and dibosons (WW, WZ, and ZZ). 

The backgrounds are subdivided into seven independent categories distinguished by their 
yields and shapes in the signal region. The six samples: tt, W+c (events with one or more 
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c jet), W+b (events with one or more b jet), W+light-flavour (u, d,s,g) jets, single top quark, 
and multijets are defined as independent categories. The small backgrounds with similar Hj 
distributions from dibosons and Z/ 7 *+jets are merged into the "Z/ y*/VV" background. Ad¬ 
ditional contributions from tt+W and tt+Z are considered negligible. All Hj distributions are 
taken from simulation. 

For the backgrounds which contribute little to the signal region (single top quark, diboson, 
Z+jets, and multijet production), the normalizations are taken directly from the simulation. 
For the four remaining processes which provide most of the background in the signal region 
(tt production, W+c, W+b, and W+light-flavour jets), the normalization is initially taken from 
simulation, but is then determined by a simultaneous fit of the background distributions to the 
data. The normalization is allowed to float freely during the limit setting. Thus, the fit finds the 
best values for these normalizations, derived using simulated and observed yields from both 
the control and signal regions. The values obtained for these normalizations for the electron 
(muon) channel are 1.01 (1.01) for tt, 2.06 (1.62) for W+c, 1.90 (1.48) for W+b, and 1.18 (1.01) 
for W+light-flavour jets. The tt background dominates and constitutes 80% of events with 1 
b-tagged jet and 93% of events with 2 or more b-tagged jets, while W+c and W+b backgrounds 
contribute to 8% and 2%, respectively. Differences in normalizations between electron and 
muon channels are accounted for in the systematic uncertainties, as noted in Table [9j 

A closure test is performed to assess the validity of the assumption that the normalizations de¬ 
rived from the fit to data are not dependent on the jet multiplicities of the samples. A sample of 
events with at least four jets, none of which are b-tagged, is used for the closure test. The agree¬ 
ment between observed and predicted events, using the post-fit values of the normalizations, 
across all bins in the high jet multiplicity region is found to be within 10%. 

8.3 Event yields 

The number of data events after different selection cuts are compared to expectations from SM 
backgrounds and are shown in Fig. [6] for both the electron and muon channels. Results are in 
good agreement with SM background expectations. 

The number of expected events in the final selection for each subsample can be seen in Ta¬ 
ble [5] The number of events for data, SM background processes, and a charged Higgs boson 
with a mass of m H + = 250 GeV are shown. The leading contributions to the SM background 
come from tt events with a semi-leptonic final state, W boson production in association with 
heavy-flavour jets, and single top quark production. Statistical and systematic uncertainties 
are evaluated as described in Section [9j 

The Hj distributions for the two signal regions in the muon channel are shown in Fig. [7j Limits 
on the production cross section of the charged Higgs boson are extracted by exploiting these 
distributions. 


9 Systematic uncertainties 


The uncertainties common to the analyses are presented in Section 9.1 The uncertainties spe¬ 


cific to the individual analyses are discussed in Sections 9.2-9.5 


9.1 Uncertainties common to the analyses 

The sources of systematic uncertainties common to the analyses (unless specified otherwise) 
and affecting simulated samples only are as follows: 
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Figure 6: Event yields after different selection cuts for both the e+jets (left) and y+jets (right) 
final state. Expectations for the charged Higgs boson for m H + = 250 GeV, for the H + —> tb 
decays, are also shown. For illustrative purposes, the signal is normalized, assuming 13 (t —> 
H + b) = 1, to a cross section of 1 pb, which is typical of the cross section sensitivity of this 
analysis. The bottom panel shows the ratio of data over the sum of the SM backgrounds with 
the total uncertainties. 


Table 5: Number of expected events for the SM backgrounds and for signal events with a 
charged Higgs boson mass of = 250 GeV in the f+jets final states after the final event 
selection. Normalizations for W+light-flavour jets, W + c, W + b, and tt are derived from data. 
Normalizations for other backgrounds are based on simulation. For illustrative purposes, the 
signal is normalized, assuming 13 (t —> H + b) = 1, to a cross section of 1 pb, which is typical of 
the cross section sensitivity of this analysis. Statistical and systematic uncertainties are shown. 


Source 

tag — 1 

tag — ^ 

Nb tag — 1 

N> tag ^ 2 

e+jets 

/(+jets 

H+ -+ tb, m H + = 250 GeV 

315 ± 4 ± 17 

647 ± 6 ± 34 

348 ± 5 ± 19 

707 ± 7 ± 37 

tt 

64111 ± 74 ± 5174 

51059 ± 66 ± 4679 

71593 ± 78 ± 5711 

57094 ± 70 ± 5160 

W+c 

8031 ± 89 ± 1047 

482 ± 21 ± 79 

7156 ±77 ±11193 

460 ± 18 ± 92 

W+b 

4470 ± 61 ± 1206 

1486 ± 35 ± 404 

3926 ± 53 ± 1386 

1364 ± 32 ± 484 

W+u,d,s,g 

3326 ± 44 ± 598 

90 ± 7 ± 21 

3231 ± 39 ± 581 

95 ± 7 ± 22 

Single top quark 

4059 ± 42 ± 463 

2253 ± 30 ± 274 

4496 ± 44 ± 524 

2493 ± 32 ± 295 

Z/7*/W 

1492 ± 54 ± 771 

237 ± 21 ± 130 

1792 ± 60 ± 942 

269 ± 22 ± 140 

Multijet background 

990 ± 270 ± 1040 

280 ± 160 ± 290 

1220 ± 480 ± 1260 

59 ± 34 ± 60 

Total SM backgrounds 

86480 ± 310 ± 5620 

55890 ± 190 ± 4720 

93410 ± 500 ± 6240 

61836 ± 95 ± 5194 

Data 

86580 

59637 

92391 

65472 
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19.7 fb' 1 (8 TeV) 



19.7 fb' 1 (8 TeV) 



Figure 7: The Hj distributions observed in data and predicted for signal and background in 
the /f+jets channel with tag = 1 (left) and Nb tag > 2 (right). Normalizations for tt, W + c, 
W + b, and W+light-flavour jets are derived from data. Normalizations for other backgrounds 
are based on simulation. Expectations for the charged Higgs boson for 7JZ h + = 250 GeV, for the 
H + —> tb decays, are also shown. For illustrative purposes, the signal is normalized, assuming 
B(t —> H + b) — 1, to a cross section of 1 pb, which is typical of the cross section sensitivity of 
this analysis. The bottom panel shows the ratio of data and the sum of the SM backgrounds 
with the total uncertainties. Bin contents are normalized to the bin width. 
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• Uncertainties in the lepton trigger, identification, and isolation efficiencies are calcu¬ 
lated from independent samples with a "tag-and-probe" method. The uncertainties 
in the single electron, single muon, and dilepton triggers amount to 2%, 2%, and 3%, 
respectively. For the r h +jets final state, the treatment is detailed in Section [972) 

• The uncertainty in the efficiency and identification of electrons is 2% (1%) for pj > 

20 (30) GeV. For muons, the uncertainty in the efficiency and identification is 1%; 

• The uncertainty in ip identification efficiency is estimated to be 6% Il83l ; 

• The misidentification uncertainty in events with an electron misidentified as the Ty 
is 20% (25%) for the barrel (endcap); for events with a muon (jet) misidentified as 
the T| 1 an uncertainty of 30% (20%) is estimated If83l : 

• The uncertainty in the ry energy scale (r h ES) is estimated by varying the ry momen¬ 
tum by ±3% j83j; 

• The uncertainties in the jet energy scale (JES), jet energy resolution (JER), and the 
contribution to Ey' ss scale from particles not clustered to jets ("unclustered Ey ,ss 
scale") are estimated independently according to the prescription described in Ref. 1421 . 
and found to within 1-6% for the signal and dominant simulated backgrounds in all 
the analyses. The variations of these quantities are also propagated to the E™ ISS . The 
uncertainty in JES is evaluated as a function of jet pj and jet //, and takes into account 
JES variations due to parton flavour; 

• The uncertainty arising from b tagging/mistagging efficiencies is estimated accord¬ 
ing to the description in Ref. Il44l . Values of 3-20% are found in the different analyses; 

• A 100% uncertainty is assumed for the reweighting of the top quark pj spectrum of 
each top quark in simulated SM tt events, discussed in Section [4} The reweighting 
and uncertainty depends on the top quark decay Ii69l ; 

• The uncertainty in pileup event modelling is estimated by varying the total inelastic 
cross section used to infer the pileup distribution in data by ±5%; 

• Uncertainties in the theoretical cross section normalization described in detail in Sec¬ 
tion |4j 

• For the /iry, E+jets, and EE' final states, the uncertainties due to ME and parton 
shower (PS) matching, and those due to the factorization and renormalization scale 
choices are applied only to the dominant simulated tt backgrounds; they are esti¬ 
mated by varying by a factor of two the threshold between jet production at the ME 
level and via PS and by varying by a factor of four the nominal scale given by the 
momentum transfer of the hard process (Q 2 ) in the event; 

• For the /Yry and EE' final states, the uncertainty in the b-tagged jet multiplicity distri¬ 
bution shapes due to PDF variations is estimated separately for the dominant simu¬ 
lated tt backgrounds by varying independently the components of the PDF param¬ 
eterization; 

• For the }IT^ and EE' final states, the uncertainty due to the modelling of the associated 
heavy-flavour production (tt+bb) is taken into account by assigning to each bin of 
the b-tagged jet multiplicity distribution of the tt+bb events an uncorrelated bin-by- 
bin uncertainty of 44%. This uncertainty is based on the comparison between the 
observed and predicted ratios of cr(tt + bb) /cr(tt + qq) HU; 

• The uncertainty in the integrated luminosity is estimated to be 2.6% Il85l . 
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9.2 The Th+jets final state for H+ —> t + v t 


In the Th+jets final state, some of the systematic uncertainties related to simulated samples also 
affect the background measurements from data. In the multijet background, a small number of 
simulated EW+tt events is subtracted from the data to obtain the number of multijet events. The 
uncertainties affecting this small number of simulated events are taken into account, but their 
magnitudes are suppressed because they apply to only a fraction of the multijet background 
and a minus sign is assigned for them to denote anticorrelation. For the "EW+tt with Th" 
background, uncertainties related to the simulated r lepton decays are taken into account. 

In addition to the uncertainties already described in Section |9.1[ the following sources of sys¬ 
tematic uncertainties are taken into account for the Th+jets final state: 


• The uncertainties in the efficiencies of the t part and E™ ISS part of the t+E™ ISS trigger 
measured from data and simulation are considered separately. The simulated sam¬ 
ples are affected by both sources of uncertainty, while the "EW+tt with Th" back¬ 
ground, obtained with the "embedding" procedure, is affected only by the uncer¬ 
tainty in the trigger efficiency measured in data. Furthermore, for the "EW+tt with 
Th" background, the data part of the p trigger efficiency is also considered, and a fur¬ 
ther 12% uncertainty is applied for approximating the Ey ISS of the high-level trigger 
by offline calorimeter-based Ey ISS ; 

• The uncertainty in vetoing events with electrons and/or muons affecting only the 
simulated samples is estimated from the uncertainty in the electron and muon recon¬ 
struction, identification, and isolation efficiencies as 2% (1%) for electrons (muons); 

• A 50% normalization uncertainty for the nij distribution is assigned for the simu¬ 
lated single top quark samples in the "EW+tt no Th" background for assigning as 
event weight the probability to pass b tagging instead of applying the b tagging 
condition; 


The uncertainties in the "EW+tt with Th" background measurement method are 
described in the following. The uncertainty in the muon identification efficiency 
in data is found to be small. The contamination of the p+jets control sample by 
multijet events is estimated with a // enriched simulated multijet sample to be at 
most 2%, which is taken as a systematic uncertainty. The fraction of events with 
W 


TVt -a \IV } 1 V T , 


discussed in Section 5.2.1 


is evaluated from simulated events 
and found to obey a functional form (1 — a) p T b , where a and b are positive con¬ 
stants and pj is the transverse momentum of the selected muon. The systematic 
uncertainty for correcting the event yield for this effect amounts to 1.2%. A 100% 
uncertainty is assumed on the event weights accounting for the difference between 
the T+jets and embedded p+jets events from simulated tt events (denoted as "Non- 
emb. vs. emb. difference" in Table [6) observed in the nij distribution; 

The uncertainties in the multijet background measurement method are described in 
the following. The statistical uncertainty in th e Ey ISS template fit that is performed 
in each bin of py h , as described in Section 


5.2 


is estimated to be 3% in each p T h bin. 
The difference in the nij distribution shapes between the nominal sample and the 
sample with inverted q, isolation criterion is taken as a systematic uncertainty. It 
is evaluated from the ratio of the event yields of the samples with nominal and in¬ 
verted ^ isolation criterion as a function of nij after requiring the other q, selection 
criteria, the veto against electrons and muons, at least three jets, and the require¬ 
ment on R™y. The statistical uncertainty of the ratio of the event yields is found 
to account for the difference in the shape and its magnitude is taken as the system- 
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atic uncertainty. Its value ranges between 5-15% depending on the bin of the nij 
distribution. 

A summary of the systematic uncertainties is shown in Table [ 6 ] 

In the region where the background yields are taken from the exponential fit on nij, the sta¬ 
tistical uncertainties in the background distributions are given by the uncertainties on the fit 
parameters while the relative values of the systematic uncertainties are kept the same like in 
the unfitted nij distribution. 

The dominant systematic uncertainties for signal arise from Th identification, Th energy scale, 
b tagging, and the theoretical tt cross section uncertainty for m H + < m t — m b- For the back¬ 
grounds, the dominant uncertainties are those in r h identification, jet—» r h misidentification, 
treatment of the £™ lss part of the trigger, and the difference between the transverse mass shapes 
of the r+jets and embedded y+jets events. In the region i > 300 GeV the sensitivity of the 
analysis is driven solely by the signal acceptance and the uncertainties in the signal. 


9.3 The ;/r h final state for H+ -a t + v t and H+ -* tb 


The dominant sources of systematic uncertainties are the Th identification and misidentifica¬ 
tion, the top quark pj modelling, and the prediction of the tt cross section. In addition to 
the uncertainties described in Section 9.1 an uncertainty associated with the misidentified Th 
background estimated from data is evaluated as half of the maximum variation between the 
"W+jet" and "multijet" estimates discussed in Section [G2] The statistical uncertainty associated 
with the number of events in the control region to which the final estimate is applied amounts 
to 1 % and is taken into account in the limit computation. 


The systematic uncertainties for the signal and background samples are summarized in Table[7] 
The diboson and Drell-Yan background yields are small compared to the uncertainty on the tt 
background, and consequently are not used in the limit computation. Results are not sensitive 
to the inclusion of those backgrounds. 


9.4 Dilepton (ee/e^/////) final states for H+ —> x + v r and H+ —* tb 

The main sources of systematic uncertainties are the unclustered £™ ISS scale, the b tagging effi¬ 
ciency, and the prediction of the tt cross section. 

The systematic uncertainties for signal and background events are summarized in Table [ 8 ] The 
diboson, Z/ 7 *, "other tt", and W+jets backgrounds yields are small compared to the uncer¬ 
tainty on the tt background, and consequently are not used in the limit computation. Results 
are not sensitive to the inclusion of those backgrounds. 


9.5 Single-lepton (e/^+jets) final states for H+ —» tb 

In addition to the uncertainties described earlier in this section, the following systematic un¬ 
certainties specific to the £+jets final states, affecting the simulated samples only, are as follows: 

• The normalizations for tt, W+c, W+b, and W+light-flavour backgrounds are left un¬ 
constrained. Statistical and systematic uncertainties are applied to yields in the con¬ 
trol regions described in Section [3] These uncertainties are based on deviations of 
the fitted normalization factor when varying multijet and Z/ 7 *+jets contributions 
by a factor of two, signal contamination by a factor of five, and by requiring either 
two or three jets in the control region. The total uncertainty in the normalization 
factors ranges between 5-35%. 
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Table 6: The systematic uncertainties (in %) on event yields for the charged Higgs boson sig¬ 
nal processes tt —> bH + bH~ (H + H~), tt —> bH + bW~ (H + W~), and pp —* t(b)H + (H + ) and 
for the background processes. The uncertainties which depend on the nij distribution bin are 
marked with (S) and for these the maximum integrated value of the negative or positive vari¬ 
ation is displayed. Empty cells indicate that an uncertainty does not affect the sample. The 
uncertainty values within the rows are considered to be fully correlated and the values within 
the columns are considered to be uncorrelated. A minus sign in front of an uncertainty value 
means anticorrelation with other values in the same row. 



Signal 

Signal 

Signal 

Multi¬ 

EW+tt 

EW+tt 


H+H- 

H+W- 

H+ 

jets 

with +h 

no T h 

r part of trigger (data) 

1.5-1.8 

1.3-1.5 

1.8-3.0 

-0.5 

1.2 

1.4 

r part of trigger (simulation) 

0.7-0 .8 

0 . 6 - 0 .7 

0 . 8 - 1.1 

- 0.2 


0.8 

E™ ss part of trigger (data) 

2.6-3.3 

25-2.8 

2.9-4.2 

- 1.2 

2.5 

2.8 

E” lss part of trigger (simulation) 
Approximation in Eip lss part of trigger 

0.1 

0.1 

0.1 

- 0.1 

12 

0.4 

Single }i trigger; data 

Veto of events with e 

0 . 1 - 0.2 

0.2-0.3 

0.2-0.3 

<- 0.1 

- 0.1 

0.4 

Veto of events with }i 

0.1 

0 . 1 - 0.2 

0.1 

<- 0.1 


0.5 

T h identification (S) 

6.0 

6.0 

5.9-6.0 

- 0.8 

6.0 


e misidentification as Th (S) 

< 0.1 

< 0.1 

< 0.1 

- 0.1 


3.3 

}i misidentification as q, (S) 

< 0.1 

< 0.1 

< 0.1 

<- 0.1 


1.1 

Jet misidentification as Tp, (S) 

0.1 

0.1-0.3 

0.1 

-6.9 


17 

energy scale (S) 

0.3-2 .6 

2.7-5.2 

0.3-2.7 

- 1.8 

5.8 

2.0 

Jet energy scale 

2.6-5.2 

2.0-3.0 

1 . 6 - 2.1 

-1.4 


3.2 

Jet energy resolution 

1 . 1 - 1.8 

0.5-1.3 

0.7-1.5 

- 0.2 


3.2 

Unclustered £™ ss energy scale 

0.1-0.4 

0.1-0.9 

0.1-0.4 

-0.5 


1.5 

b-jet tagging (S) 

5.9-20 

4.7-5.3 

4.6-54 

-3.5 


5.0 

Top quark pj modelling (S) 




+5.6 

-6.8 


+n 

—6.6 

Pileup modelling 
}i identification; data 

0.1-0.9 

0 . 1 - 0.8 

0 . 1 - 0.6 

- 0.1 

<- 0.1 

2.9 

Multijet contamination 





2.0 


W — > tv t — > pv^v-r fraction 





1.2 


Non-emb. vs. emb. difference (S) 
Multijet nij distribution shape (S) 




4.6 

+14 

-12 


Multijet template fit 

Probabilistic nij in single top quark 




3.0 


6.8 

tt cross section, scale 

+2.5 

-3.4 

+2.5 

-3.4 


+1.0 

-0.7 


+2.2 

-2.9 

tt cross section, PDF+as 

4.6 

4.6 


- 1.6 


4.0 

Single top quark cross section 






1.0 

W+jets, Z/ 7 *, W cross section 






0.1 

Integrated luminosity 

2.6 

2.6 

2.6 

- 0.8 


2.6 
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Table 7: The systematic uncertainties (in %) for the }ir^ final state for backgrounds, and for 
signal events from H + —» tb decays for »z H = 250 GeV. These systematic uncertainties are 
given as the input to the exclusion limit calculation. The uncertainties that depend on the 
b-tagged jets multiplicity distribution bin are marked with (S) and for these the maximum 
integrated value of the negative or positive variation is displayed. Empty cells indicate that an 
uncertainty does not affect the sample. The uncertainty values within the rows are considered 
to be fully correlated and the values within the columns are considered to be uncorrelated. The 
uncertainties in the cross sections are to be considered uncorrelated for different samples and 
fully correlated for different final states of the same sample (e.g. the different ft decays). 



Signal 

tt —y fiTfr + X 

tt dilepton 

r h mis-id single top quark 

Single }i trigger 

2.0 

2.0 

2.0 


e identification 

2.0 

2.0 

2.0 

2.0 

}i identification 

1.0 

1.0 

1.0 

1.0 

Th identification 
e misidentification as Tj, 

6.0 

6.0 

3.0 

6.0 

}i misidentification as 



3.0 


Jet misidentification as ip 



20 


r h energy scale (S) 

0.6 

2.4 

4.4 

4.1 

Jet energy scale (S) 

2.5 

1.9 

2.6 

3.9 

Jet energy resolution (S) 

0.8 

0.1 

1.6 

0.2 

Unclustered £™ ss energy scale (S) 

0.8 

0.1 

1.8 

0.2 

b tagging (S) 

1.8 

1.8 

2.7 

3.2 

udsg—>b mistagging (S) 

<0.1 

<0.1 

<0.1 

0.1 

Top quark pj modelling (S) 


5.4 

5.2 


Pileup modelling 

4.0 

2.0 

8.0 

2.0 

Misidentified ip background 




11 

Cross sections 


+ 2 - 5 ±46 

-3.4 ^ ^- D 

+ 2 - 5 ±46 

-3.4 ^ 

8.0 

Matching scale (S) 


12 

5.1 


Fact./renorm, scale (S) 


3.4 

7.5 


PDF effect on shape 


shape only 

shape only 


Heavy flavours (S) 


<0.1 

<0.1 


Integrated luminosity 

2.6 

2.6 

2.6 

2.6 
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Table 8 : The systematic uncertainties (in %) for backgrounds, and for signal events from 
H + —> tb decays for the dilepton channels for a charged Higgs boson mass in H( = 250 GeV. The 
e/< final state is shown as a representative example. These systematic uncertainties are given as 
the input to the exclusion limit calculation. The uncertainties that depend on the b-tagged jets 
multiplicity distribution bin are marked with (S) and for these the maximum integrated value 
of the negative or positive variation is displayed. Empty cells indicate that an uncertainty does 
not affect the sample. The uncertainty values within the rows are considered to be fully corre¬ 
lated and the values within the columns are considered to be uncorrelated. The uncertainties in 
the cross sections are to be considered uncorrelated for different samples and fully correlated 
for different final states of the same sample (e.g. the different tt decay channels). 



Signal 

tt dilepton 

z/7* -> a 

single top quark 

e// trigger 

3.0 

3.0 

3.0 

3.0 

e identification 

2.0 

2.0 

2.0 

2.0 

j.i identification 

1.0 

1.0 

1.0 

1.0 

Jet energy scale (S) 

1.4 

1.1 

1.7 

1.4 

Jet energy resolution (S) 

0.3 

0.3 

0.4 

0.4 

Unclustered £™ ss energy scale (S) 

1.3 

2.1 

11.7 

2.6 

b tagging (S) 

2.4 

3.7 

10 

4.3 

udsg—S>b mistagging (S) 

2.3 

3.6 

10 

4.4 

Top quark pj modelling (S) 


3.8 



Pileup modelling 

0.6 

0.4 

1.2 

1.2 

Cross sections 


^4 ±4.6 

4.0 

8.0 

Matching scale (S) 


7.7 



Fact./renorm, scale (S) 


8.4 



PDF shape 


shape only 



Fleavy flavours (S) 


<0.1 



Integrated luminosity 

2.6 

2.6 

2.6 

2.6 



27 


• A 50% uncertainty Ii86fj88l is applied to the Z/ 7 *+jets and diboson backgrounds due 
to their small contribution to the signal region; 

• A 100% systematic uncertainty is applied to the QCD cross section normalization. 
This accounts for the maximal variation in the QCD normalization when left uncon¬ 
strained in the background-only fit to data while constraining normalizations for 
other backgrounds to their systematic uncertainties. 

The systematic uncertainties for signal and background events are summarized in Table [9] 


Table 9: The systematic uncertainties (in %) for backgrounds, and for signal events from 
H + —> tb decays for the f’+jets channels for a charged Higgs boson mass 7)z H + = 250 GeV. The 
uncertainties that depend on the shape of the Hj distribution bin are marked with (S) and for 
these the maximum integrated value of the negative or positive variation is displayed. Empty 
cells indicate that an uncertainty does not affect the sample. The uncertainty values within the 
rows are considered to be fully correlated, with the exception of cross section and data-driven 
normalization, which are considered to be uncorrelated. The uncertainty values within the 
columns are considered to be uncorrelated. Uncertainties labelled with a are only present 
in the CR with an implicit unconstrained parameter correlated across all bins (Sec. 8.21. The 
values for these are assigned prior to the setting of limits. 



H+ -> tb 

tt 

W+c 

W+b 

W+u,d,s,g 

single top quark 

Z/7VVV 

Multijets 

Single-e trigger 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

2.0 

Single-// trigger 

2.0 

2.0 

2.0 

2.0 
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10 Results 

A statistical analysis of the mj (Fig. [3]), b-tagged jet multiplicity (Fig. [4] (right) and Fig. [5] (right)), 
and Hj (Fig. [7jl distributions has been performed using a binned maximum likelihood fit. The 
data agree with the SM prediction and consequently 95% CF upper limits on charged Higgs 
boson production are derived using the modified frequentist CF S criterion Il89l |90l with a test 
statistic based on the profile likelihood ratio with asymptotic approximation [91, 92J. 

The systematic uncertainties described in Section [9] are incorporated via nuisance parameters 
following the frequentist paradigm. Correlations between the different sources of systematic 
uncertainty are taken into account. Uncertainties affecting the shape of the mj, b-tagged jet 
multiplicity, or Hj distributions are represented by nuisance parameters whose variation re¬ 
sults in a continuous perturbation of the distribution l93l . 
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10.1 Model-independent limits on charged Higgs boson production (H+ —>■ 

r+v T ) 

In the analysis of the H + —> r + v T decay mode with the ry+jets final state no assumption on the 
charged Higgs boson branching fractions is needed because subtracting the background from 
"EW+tt with iy" will remove any potential H + —> tb and other such signals from data due 
to the embedding technique described in Section 5.2.1 For m H + = 80-160 GeV, the charged 


Higgs boson is produced most copiously through tt production which can produce one (tt —> 
bH + bW~) or two charged Higgs bosons (tt —> bH + bH~) if B(t — > H + b) > 0. Furthermore, the 
presence of the charged Higgs boson suppresses the tt —> bW bW yield compared to the SM 
prediction. Consequently, the number of events in a given bin of the mj distribution depends 
on the signal strength parameter }i according to: 


N(n) = fi 2 s(H+H“) + 2/i(l - //) s(H+W _ ) + (1 - }i) 2 b(W+W“) + b, (7) 

where }i = B{ t —>• H + b) B( H + — > t + v t ), s(H + H~) and s(H + W~) are the number of expected 
signal events for the tt —> bH + bH _ and tt —» bH + bW~ processes, respectively; & (W ! W ) is 
the expected number of events from the portion of tt —> bW 1 bW background that is estimated 
with simulation, and b is the expected number of other background events. The number of 
signal and tt —> bW 1 bW background events is normalized to the SM predicted cross section 
and by setting B(t — > H + b) B (H + —> t + v t ) = 1 for a top quark decaying to a charged Higgs 
boson. 

For m H + = 180-600 GeV, the number of events in a given bin of the nij distribution depends 
on the signal strength parameter according to: 


N(}i) = }L£ S C + b, 


( 8 ) 


where }i = cr( pp —> t(b)H + ) B(H 1 —> t ' v T ), e s is the event selection efficiency for signal 
events, C is the integrated luminosity, and b is the expected number of background events. 

The upper limits on B (t —> H + b) B( H + —> t + v t ) and on cr( pp —» t(b)H + ) B( H + —> t + v t ) 

are shown in Fig. [8] for the H + —» r + v T decay mode with the ry+jets final state for the ranges 
/n H + = 80-160 GeV and m H + = 180-600 GeV, respectively. The numerical values of the limits 


are given in Table 10 At m H = 250 GeV an excess of data is observed with a local p-value of 
0.046 corresponding to significance of 1.7a. 

10.2 Limits on charged Higgs boson production with branching fraction as¬ 
sumed 

In the presence of a charged Higgs boson and for m H = 180-600 GeV, the analyses of the //Ty, 
f+jets, and IB final states have sensitivity to both H + —> t + v t and H + —> tb decays. Con¬ 
sequently, a model-independent limit can neither be provided for cr(pp —> t(b)H + ) B( H + —> 
t + v t ) nor for <r(pp —> t(b)H + ) 6(H + —> tb). Nevertheless, one can test models by fixing 
B( H + —> t + v t ) and B( H + —> tb). In this section, results are reported for a model with 
B( H + —> tb) = 1, to which the ry+jets anal ysis is blind because of the estimates of the back¬ 
grounds from data like described in Section 5.2.1 For B (H + —> r ' v r ) = 1, the sensitivity of 
the //ly and ££' final states analyses is found to be substantially weaker than that obtained in 
the Th+jets analysis. 






10.2 Limits on charged Higgs boson production with branching fraction assumed 


29 


19.7 m' 1 (8 TeV ) 



Figure 8: Expected and observed 95% CL 
H + b) £?(H+ —* t + v t ) with m H + = 80-160 GeV 
with 7 /z h + = 180-600 GeV (right) for the H + 
regions above the solid lines are excluded. 



model-independent upper limits on £>(t —> 
(left), and on c(pp —> t(b)H+) £>(H + —* t + v t ) 
t 1 v T search in the r h +jets final state. The 


Table 10: Expected and observed 95% CL model-independent upper limits on B(t —> 
H + b) £>(H+ —* r+v T ) for m H + = 80-160 GeV (top), and on rr(pp —> t(b)H+) B( H + —* z + v T ) 
for m H . = 180-600 GeV (bottom), for the H + —> r ' v r search in the r h +jets final state. 


m H+ 

[GeV] 

-2a 

Expected limit 
— la median +1(7 

+2(7 

Observed 

limit 

95% CL upper limit 

on B (t — 

► H+b) B( H+ -+ r 

f t/ T ) 

80 

0.0059 

0.0079 

0.0112 

0.0160 

0.0221 

0.0120 

90 

0.0042 

0.0057 

0.0080 

0.0115 

0.0160 

0.0092 

100 

0.0033 

0.0044 

0.0062 

0.0089 

0.0124 

0.0061 

120 

0.0018 

0.0024 

0.0034 

0.0049 

0.0069 

0.0028 

140 

0.0012 

0.0017 

0.0024 

0.0034 

0.0048 

0.0017 

150 

0.0011 

0.0015 

0.0021 

0.0031 

0.0043 

0.0015 

155 

0.0012 

0.0016 

0.0023 

0.0033 

0.0046 

0.0016 

160 

0.0011 

0.0016 

0.0022 

0.0032 

0.0045 

0.0015 

95% CL upper limit on cr (pp —> t(b)H+) B( H + —)■ r + v T ) [pb] 

180 

0.213 

0.289 

0.409 

0.587 

0.816 

0.377 

190 

0.188 

0.254 

0.358 

0.516 

0.719 

0.373 

200 

0.152 

0.205 

0.291 

0.423 

0.587 

0.361 

220 

0.114 

0.155 

0.221 

0.321 

0.448 

0.332 

250 

0.081 

0.110 

0.159 

0.231 

0.328 

0.267 

300 

0.048 

0.065 

0.096 

0.142 

0.205 

0.153 

400 

0.022 

0.032 

0.049 

0.076 

0.115 

0.054 

500 

0.014 

0.021 

0.033 

0.056 

0.088 

0.032 

600 

0.011 

0.016 

0.028 

0.047 

0.076 

0.025 
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Equation (|8j is used to derive the limits by counting the number of events in bins of the b- 
tagged jet multiplicity distribution for the and ££' final states, and in bins of the Hj distribu¬ 
tion for the f+jets final state. The upper limits on <r(pp —> t(b)H + ) assuming B( H + —> tb) = 1 
are shown in Fig. |9]for the //r h (top left), f'+jets (top right), and ££! (bottom) final states. 


19.7 fb" 1 (8 TeV) 



nv [GeV] 



"V [GeV] 



m H * [GeV] 


Figure 9: Expected and observed 95% CL upper limits on u(pp —> t(b)H + ) for the /rq, (upper 
left), f+jets (upper right), and ££' final states (bottom) assuming B (H ' —>■ tb) = 1. The regions 
above the solid lines are excluded. 


The upper limit on <r(pp 


t(b)H + ) for the combination of the «T h , f+jets, and ££' final states 

In the combination, the 


is shown in Fig. 10 The numerical values are reported in Table 11 


sensitivity is driven by the £+jets final state. 


10.3 Combined limits on tan/3 in MSSM benchmark scenarios 

Using all decay modes and final states, exclusion regions have been set in the »z H +-tan (1 plane 
according to the LHC Higgs cross section working group prescription for different MSSM 
benchmark scenarios l|29l |33|: "updated m™ ax ", "m™ od+ ", "m™ od '", "light stop", "light stau", 
"tau-phobic", and "1ow-Mh" scenarios. These MSSM benchmark scenarios are compatible 
with the properties of the recently discovered neutral scalar boson and with the current bounds 
on supersymmetric particle masses, and they are specified using low-energy MSSM parame¬ 
ters, i.e. no particular soft SUSY-breaking scenario is assumed. The updated »z™ ax scenario and 
m !(' od scenarios allow the discovered scalar boson to be interpreted as the light CP-even Higgs 













10.3 Combined limits on tan/5 in MSSM benchmark scenarios 
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Figure 10: Expected and observed 95% CL upper limits on a (pp —> t(b)H + ) for the combina¬ 
tion of the }i Th, f+jets, and U final states assuming B{ H + —» tb) = 1. The region above the 
solid line is excluded. 


Table 11: Expected and observed 95% CL upper limits on tr(pp —> t(b)H + ) B( H + —> tb) as¬ 
suming B{ H + —5 tb) = 1 for the combination of the /( Th , f+jets, and lit' final states. 

m H + Expected limit [pb] Observed limit [pb] 

[GeV] —2m —In median +1(7 +2(7 limit 


95% CL upper limit on a (pp —> t(b)H + ) with B( H + -+ tb) = 1 


180 

1.07 

1.43 

2.01 

2.81 

3.78 

1.99 
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0.87 

1.16 

1.62 

2.27 

3.07 

1.52 
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0.62 

0.83 
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1.64 
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0.99 
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0.49 

0.66 

0.93 

1.31 

1.78 

0.89 
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0.33 

0.45 

0.62 

0.88 

1.18 

0.54 

400 

0.22 

0.29 

0.40 

0.57 

0.76 

0.33 

500 

0.15 

0.20 

0.28 

0.39 

0.52 

0.21 
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0.10 

0.14 

0.19 

0.27 

0.36 

0.13 
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boson in large parts of the m H i -tan p plane. The light stop scenario leads to a suppressed rate 
for the Higgs boson production by gluon fusion, and the light stau scenario enhances the decay 
rate of the light CP-even Higgs boson to photons. A tau-phobic scenario has suppressed cou¬ 
plings to down-type fermions. In the 1ow-Mh scenario, the discovered scalar boson is assumed 
to be the heavy CP-even Higgs boson and 7 «a is fixed to be 110 GeV causing m H + to be 132 GeV. 


Figure 11 shows the limits on the updated m™ ax and m j) lod " scenarios. For m H + = 90-160 GeV, 
the analysis of the H + —* t + v t decay mode with the +jets final state described in Section[5]is 
taken as input. The mass range starts here from m H + = 90 GeV, as the lower values of a charged 
Higgs boson mass are not accessible in the considered MSSM scenarios. For m H + = 200- 
600 GeV, a combination of all decay modes and final states is used to set the limits. In this com¬ 
bination, the signal yields from the H + —> r + v T and H + —> tb decay modes are defined by the 
branching fractions predicted by the model. If the limit on the charged Higgs boson production 
for a given m H -tan /3 point is smaller than the cross section predicted by the model II281 - 13H . 
the point is excluded. The mass range is chosen to start from m H + = 200 GeV to avoid the 
interference region where a charged Higgs boson is produced both from off-shell top quark 
decays and through direct production. In all these scenarios except for the 1 ow-Mh and light 
stop scenarios, a lower bound of about 155 GeV on the charged Higgs boson mass has been set 
assuming m ^ = 125 ± 3 GeV. The light stop scenario is excluded for m H + < 160 GeV assuming 
mh = 125 ± 3 GeV. For m H + > m t — rn b , the H + —> tb decay mode searches yield a lower 
limit on tan (5 while the upper limit on tan j 6 is dominated by the results from the analysis of 
the H + —> r 1 i' T decay mode with the ih+jets final state. The 1 ow-Mh scenario is completely 
excluded (Fig. 121 assuming the heavy CP-even MSSM Higgs boson mass is rn\\ = 125 ± 3 GeV. 


In Figs. 11 -[12} theoretical systematic uncertainties affecting the expected signal event yields 
are added to the limit computation, modelled as nuisance parameters, in addition to the un¬ 
certainties discussed in Section [9] The uncertainty in the branching fractions of the charged 
Higgs boson is estimated from the decay width uncertainties as in Ref. 11941 by scaling each 
partial width separately while fixing all others to their central values. This results in individ¬ 
ual theoretical uncertainties for each branching fraction. The width uncertainties comprise the 
uncertainty from missing higher order corrections to beyond LO EW diagrams (5%), missing 
higher order corrections to NLO QCD (2%), and -correction uncertainties (3%) Il95l . The A 
correction arises from the presence of squarks and gluino contributions in the charged Higgs 
boson Yukawa coupling to top and bottom quarks 


For m H i = 90-160 GeV, the theoretical uncertainties in the signal yield include the uncertain¬ 
ties in the branching fractions for t —» H + b and H + —» r 1 v T totalling 0.1-5.0% depending on 
m H + and tan /I. Additionally, an uncertainty of 3% is added to the simulated tt background to 
take into account higher order corrections to the tt cross section. For m H + = 200-600 GeV, the 
charged Higgs boson production cross section uncertainty and the uncertainty in the branch¬ 
ing ratios are considered. The cross section uncertainty varies between 22-32% depending on 
m H +, tan (5, and the MSSM benchmark scenario. The uncertainty in £>(H + —> t + i/ t ) varies 
between 0.4-5.0% for tan p = 10-60 depending on 77Z H + and the MSSM benchmark scenario. 
The £>(H + —> tb) uncertainty varies between 0.1-5.0% for tan p = 1-10 depending on 77Z H + and 
the MSSM benchmark scenario. The theoretical branching fraction uncertainties for a given 
7tz H +-tan j 6 point are summed linearly according to the FHC Higgs cross section working group 
prescription [94, 95j, but the cross section and branching fraction uncertainties are treated as 
independent nuisances. The expected limit improves by no more than 2% if the theoretical 
uncertainties are treated in the statistical model as independent sources. 
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Figure 11: Exclusion region in the MSSM m H +-tan/3 parameter space for m H + = 80-160GeV 
(left column) and for ??z H + = 180-600 GeV (right column) in the updated MSSM m™ ax scenario 
(top row) and m|( lod scenarios I29l[33j (bottom row). In the upper row plots the limit is derived 
from the H ! —> r + v T search with the T^+jets final state, and in the lower row plots the limit 
is derived from a combination of all the charged Higgs boson decay modes and final states 
considered. The iluand ±2<r bands around the expected limit are also shown. The light-grey 
region is excluded. The red lines depict the allowed parameter space for the assumption that 
the discovered scalar boson is the lightest CP-even MSSM Higgs boson with a mass = 125 ± 
3 GeV, where the uncertainty is the theoretical uncertainty in the Higgs boson mass calculation. 


11 Summary 

A search is performed for a charged Higgs boson with the CMS detector using a data sam¬ 
ple corresponding to an integrated luminosity of 19.7 ± 0.5 fb” 1 in proton-proton collisions at 
s = 8 TeV. The charged Higgs boson production in tt decays and in pp —> t(b)H + is studied 
assuming H + —> r + v T and H + —> tb decay modes, using the ih+jets, }i Th, f+jets, and K:' final 
states. Data are found to agree with the SM expectations. 

Model-independent limits without an assumption on the charged Higgs boson branching frac¬ 
tions are derived for the H + —» r 1 v T decay mode in the Th+jets final state. Upper limits at 95% 
CL of B {t -A H+b) B{ H+ ^ r+v T ) = 1.2-0.15% and n(pp t(b)H+) B( H+ ^ t+v t ) = 0.38- 
0.025 pb are set for charged Higgs boson mass ranges m H + = 80-160 GeV and m H i = 180- 
600 GeV, respectively. 
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Figure 12: Exclusion region in the MSSM Higgsino mass parameter (/<) vs. tan /3 parameter 
space in the 1ow-Mh scenario |29l 1331 with rn\ = 110 GeV for the H + —> r 1 v T search with 
the Th+jets final state. The light-grey region is excluded and the blue region is theoretically 
inaccessible. The area inside the red lines is the allowed parameter space for the assump¬ 
tion that the discovered scalar boson is the heavy CP-even MSSM Higgs boson with a mass 
= 125 ± 3 GeV, where the uncertainty is the theoretical uncertainty in the Higgs boson 
mass calculation. 


Assuming £>(H + —* tb) = 1, a 95% CL upper limit of n(pp —* t(b)H + ) = 2.0-0.13pb is set for 
a combination of the jiT^, 1+jets, and W final states for m H + = 180-600 GeV. This is the first 
experimental result on the H + —> tb decay mode. Here, cross section <j(pp —> t(b)H ± ) stands 
for the sum <r(pp —> t(b)H + ) + cr( pp —> t(b)H~). 

The results are interpreted in different MSSM benchmark scenarios and used to set exclusion 
limits in the m H +-tan /I parameter spaces. In the various models, a lower bound on the charged 
Higgs boson mass of about 155 GeV is set assuming = 125 ± 3 GeV. The light-stop scenario 
is excluded for m H + < 160 GeV assuming = 125 ± 3 GeV, and the 1ow-Mh scenario defined 
in Refs. [29 331 is completely excluded assuming /«h = 125 ± 3 GeV. 
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